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We analyze atom-atom correlations in the s-wave scattering halo of two colliding condensates. 
By developing a simple perturbative approach, we obtain explicit analytic results for the collinear 
(CL) and back-to-back (BB) correlations corresponding to realistic density profiles of the colliding 
condensates with interactions. The results in the short time limit are in agreement with the first- 
principles simulations using the positive-P representation and provide analytic insights into the 
experimental observations of Perrin et al. [Phys. Rev. Lett. 99, 150405 (2007)]. For long collision 
durations, we predict that the BB correlation becomes broader than the CL correlation. 
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Experiments with colliding Bose-Einstein condensates 
(BECs) |l|, |2fl are currently attracting considerable at- 
tention in the field of ultracold quantum gases [3, U, la, 
Is III la Isi lllA 111!- A recent breakthrough in this area 
is the direct detection |2j of atom-atom pair correlations 
in the s-wave scattering halo formed in the collision of 
metastable helium ( 4 He*) condensates. Experimental ad- 
vances like this pose increasingly demanding challenges 
to theory due to the need to provide quantitatively ac- 
curate descriptions in realistic parameter regimes. 

Theoretical developments are taking place on two 

fronts using (i) using numerical techniques based on 

stochastic phase-space methods (such as the first- 

rinciples simulations in the positive-P representation 



ll| or approximate simulations based on the truncated 



Wigner- function expansion |8(), and (ii) analytic meth- 
ods such as the Bogoliubov theory |a, la |7| , undepleted 
source approximation [j, |ll|, and a simple Gaussian 
ansatz [10|. (The analytic methods often represent ex- 
tensions of quantum optics approaches describing closely 
related systems of parametric down-conversion and four- 
wave mixing 12].) Despite these developments, the exist- 
ing numerical techniques still fall short of fully describing 
the experimental measurements of Ref. |2J, whereas the 
approximations of the analytic methods are usually too 
severe to lead to full quantitative agreement with the ex- 
periments. 

In this paper, we develop an alternative analytic 
scheme to obtain explicit results for atom-atom corre- 
lations in condensate collisions. The scheme is rather 
simple, yet it compares surprisingly well with the exact 
positive-P simulations in the short-time limit. The main 
advantages of the method are the analytic transparency 
and the fact that it can model realistic density profiles 
of the colliding condensates with interactions. This is 
important for addressing the role of mode mixing due 
to the inhomogeneity of trapped atomic clouds, which is 
the strongest effect that influences the strengths and the 
width of atom-atom correlations. Our results may also 
have implications for future experiments aimed at pro- 



ducing relative atom number squeezing 11. 13 J . Jl4L Il5j ] 
and Einstcin-Podolsky- Rosen correlations 16l.ll7|. 

Additionally, we perform first-principles positive-P 
simulations of the collision dynamics and analyze the 
widths of the collinear (CL) and back-to-back (BB) corre- 
lations as per measurements of Ref. |2] . These character- 
ize, respectively, the pair correlations between the scat- 
tered atoms propagating in the same and in the opposite 
directions. The simulations are performed for 23 Na atoms 
(instead of 4 He* used in |2[) as this case appears to have 
more favorable parameters for the positive-P method to 
remain valid for long collision durations. For compari- 
son, in the case of 4 He* [2, LUj , the simulation durations 
were much shorter than the experimental collision time 
due the smaller mass and larger scattering length of 4 Hc* 
atoms. The surprising result that we find here is that the 
width of the BB correlation grows with time and eventu- 
ally becomes larger than the width of the CL correlation. 
This observation is in agreement with the measured cor- 
relation widths of Ref. 11]; it is also accompanied by 



the reduction of the BB correlation strength below the 
CL correlation, which in turn implies absence of relative 
number squeezing in the long-time limit. 

We start by considering the equations of motion de- 
scribing the collision of two BECs in the Bogoliubov ap- 
proximation in which the atomic field operator is split 
into the mean- field (MF) and fluctuating components, 
*(x,t) = *+Q(x,t)+*_o(x,t)+?(x,t). Here, *±g(x,t) 
represent the MF amplitudes of the colliding condensates 
(which can be created by splitting a single stationary con- 
densate using a Bragg pulse or Raman lasers) with mean 
momenta +Q and — Q along the x axis (in wave num- 
ber units, in the center-of-mass frame), whereas <5(x, t) 
is the fluctuating component, which is treated quantum 
mechanically to the lowest order in perturbation theory. 
The fluctuating component satisfies the equation of mo- 
tion 



ih 



d( 5(x, t) 
dt 



2m 



?(x, t)+2U%Q(x, *)«_q(x, 0(x, t), 

(1) 



where U = 4irh 2 a/m is the coupling constant describing 
s-wave scattering interactions, with a being the scatter- 
ing length. The MF components, on the other hand, sat- 
isfy the standard time-dependent Gross-Pitaevskii (GP) 
equation. However, if the mean kinetic energy of the 
colliding atoms is much larger than the interaction en- 
ergy and if we restrict ourselves to short collision times, 
then the MF amplitudes can be approximated by freely 
propagating fields [6|: v 1'±q(x, t) = y / p(x)/2 exp(±iQx — 
ihQ 2 t/2m). Here, p(x) is the density profile of the ini- 
tial source condensate and we have additionally assumed 
that the center-of-mass displacements x =F v^t (where 
v = hQ/m is the collision velocity) of the MF com- 
ponents are negligible in the short-time limit. This is 
justified if the displacements during the effective colli- 
sion time are much smaller than the characteristic size 
of the colliding BECs in the collision direction. In the 
experiments of Ref. [2| this was indeed the case. 

Substituting the expressions for \I>±q(x, t) into Eq. (JTJ 
and transforming to a rotating frame S(x, t) — v 
(5(x, t) exp(— ihQ 2 t/2m), the equation of motion for 
<5(x, t) can be rewritten in a simpler form, 



dS{x.,t) 
dt 



w 2 

2m 



HQ 2 



2m 



?(x,i)-t S (x)?(x,i), (2) 



where <?(x) is a spatially dependent effective coupling 
given by <?(x) = 2f/* +Q (x, 0)*_q(x, 0)/h = Up(x)/H. 

Converting to Fourier space, <5(x, t) = 

J dka(k, t) exp(ik • x)/(27r) 3 / 2 , yields the following 
equation of motion for the amplitude operator a(k, t) : 



dd(k, t) 
dt 



-zAfca(k, t) — i 



rfq 



(2tt) 3 /- 



r g(q + k)at( q ,i). (3) 



Here g(k) = J dxe lk ' x <7(x)/(27r) 3 / 2 is the Fourier trans- 
form of g(x), A fe = h(k 2 - Q 2 )/ (2m), and k 2 = |k| 2 . 

From Eq. ([3|) we can easily recognize the role of mode 
mixing in the spatially inhomogeneous treatment com- 
pared to an idealized uniform treatment. In the present 
inhomogeneous case, the finite width of the coupling g(k) 
implies that a(k) couples not only to a'(— k), but also to 
a range of momenta around — k, within — k ± <5k. The 
spread in <5k determines the width of atom-atom corre- 
lations with nearly opposite momenta and is ultimately 
related to the momentum width of the colliding BECs. 

To quantify the pair correlations expected between the 
s-wave scattered atoms with equal but opposite momenta 
due to momentum conservation and between the atoms in 
the collinear direction due to quantum statistical effects 
we use Glauber's second-order correlation function 

(2) _ (qt(k,t)at(k / ,t)a(k',t)a(M)) 

9 ^ k ^>- n(k,t)n(k',t) ' W 

which describes the density-density correlations between 
two momentum components k and k'. The normalization 



with respect to the product of the densities n(k, t) and 
n(k',i) ensures that </ 2 ^(k, k',t) = 1 for uncorrelated 
states. The averaging is with respect to the vacuum ini- 
tial state for the scattered modes. 

The linearity of Eq. ([3]) ensures that one can apply 
Wick's theorem to Eq. (01 and factorize the fourth-order 
moment, yielding 

9 < 2)(k , k , t) = 1+ K^|l«^. (5) 

Here, n(k,k',t) = (a f (k, t)a(k', t)) and m(k,k',t) = 
(a(k, t)a(k' , t)) are the normal and anomalous densities; 
n (k, t) = n (k, k, t) is the momentum distribution. 

It follows from Eq. ^ and the shape of g(k) that, 
for sufficiently large collision momentum Q (much larger 
than the momentum spread of the colliding BECs), the 
normal density n(k, k',t) acquires nonzero population 
primarily for pairs of nearby momenta, k' ~ k, while 
the anomalous density m(k, k', t) does so for pairs of mo- 
menta that are nearly opposite, k' ~ — k. Accordingly 
we can concentrate on the CL and BB correlations, which 
are denoted via g^ (k , Ak , t) — g^(k, k + Ak, t) and 
g4e(k, Ak, t) = g (2) (k, -k + Ak, t) and are given by 



(2) 



5 ^(k,Ak,i) = I 
5 ^(k,Ak,i) = H 



|n(k,k+Ak,t)| 2 
n(k,t)n(k+Ak,t)' 

|m(k,-k+Ak,t)| 2 
n(k,i)n(-k+Ak,i)' 



(6) 

(7) 



To calculate these correlation functions in the short 
time limit, we proceed with the Taylor expansion in time, 
up to the terms of order t 2 [l5( : 



a(k,t) =o(k,0) 



da(k, t) 



dt 



t + 



d 2 a{k, t) 



8t 2 



t A 

~2 + 



The expansion is valid for t <Si to, where to = l/g(0) = 
l/[Up(0)/h] is the time scale [18|. Using the right-hand 
side of Eq. ([3]), this gives, up to the lowest-order terms, 



n(k,k',i) 



|m(k,k',i)| 



t 2 J dqg{c l + k)g(c i + k')/{2^f 

t 2 f dxe-*( k - k ')- x [.g(x)] 2 /(27r) 3 , (8) 
i|5(k + k')|/(27r) 3 / 2 



= t 



|dxe-*( k+k '>- x 5 (x)/(2^ 



(9) 



From these results we see that the width of the CL corre- 
lation, Eq. (JBJ), is determined by the square of the Fourier 
transform of the square of the effective coupling <?(x). 
The width of the BB correlation, Eq. ([7]), on the other 
hand, is determined by the square of the Fourier trans- 
form of <7(x). Therefore, the CL correlation is generally 
broader than the BB correlation at short times. 



Thomas-Fermi (TF) "parabolic density profile. - - We 
now give explicit analytic results for the case when the 
initial condensate density profile is given by the ground 
state of the GP equation in a harmonic trap in the TF 
regime: p(x) = p (l - E.zf/flf) for E^f/^l < 1 
(i = x, y, z) and p(x) = elsewhere, with Ri being 
the TF radius along direction i and po the peak density. 
For definitcncss, we consider CL and BB correlations for 
which the displacement Ak is along one of the Carte- 
sian coordinates, i.e., k' = ik + e^Afci, where ei is the 
unit vector in the ki direction. The integrals in Eqs. ([5]) 
and ([9]) can be performed explicitly, in terms of Bessel 
functions J v (z) [lj|, yielding 



?i(k, k+e^A/ci, t) 



|m(k,-k+e 2 Afci,i)|: 



8t 2 g 2 (0)R J V2 (AkjRi) 

(2 7 r)3/2(Afc lJ R l ) 7/2 : 
2tg(0)R 3 J 5/2 (Ak l R l ) 
(27r) 3 / 2 (A/c l i? l ) 5/2 ' 



(10) 
(11) 



where g(0) = Up /h, and R = (R x R y R z ) 1 / 3 is the 
geometric mean TF radius. Applying these results to 
Afc 4 = 0, using J u (z) ~ {z/2) v fY{v + 1) for z < 1 
(y 7^ — 1, —2, . . .), we obtain the following results for 
the atomic momentum distribution and the anomalous 
density: n(k,t) ~ 4t 2 g 2 (0)R 3 / I05ir 2 and |m(k, -k,t)| ~ 

tg(0)R /I5ir 2 . These are uniform in the short-time limit, 
corresponding to spontaneous initiation of the scatte ring 
which populates the scattering modes uniformly [6|, [ll( 
without the need to strictly conserve energy. A narrow 
scattering shell around |k| = Q forms later in time, while 
the momentum cutoff, which must be assumed for a 5- 
function interaction potential, prevents the total atom 
number from diverging. 

Substituting Eqs. flTUl) and (JTTJ) into Eqs. © and (0), 
and suppressing the uniform dependency on k, we obtain 
the following explicit results for the CL and BB correla- 
tions, valid for t <^to'- 



9 { cl(^k t ,t) 



ffg>(Afc,t)~l 



105 2 7T 
105 2 7T 

32£V(0) 



JyzjAkjRj) 
(AkiRiY/ 2 

J5/2{Ak l R i 



(AkiRi)^ 



(12) 
(13) 



The CL correlation shows the Hanbury Brown and 
Twiss (HBT) bunching with the peak value of <?(x(0, t) = 
2. The BB correlation, on the other hand, shows super- 

(21 

bunching (g B g(0, £) > 2) between atom pairs with equal 

but opposite momenta, with the peak value g^(0,t) = 
1 + 7 2 /[2H 2 g 2 (0)} > 1 partly reflecting the fact that typ- 
ical mode occupation numbers in this spontaneous scat- 
tering regime are much smaller than 1. From the above 
results we also determine the widths of the CL and BB 
correlations, which are defined for simplicity as the half 
width at half maximum, w\ ~ 1.23w- and w\ ~ 



1.08w;^ , giving the ratio of w\ /w\ ~ 1.14. Here 
w\ ~ 1.99/Ri is the width of the momentum distri- 
bution p(ki) = | / dx v //o(x) exp(-ifc 4 Xi)/(27r) 3 / 2 | 2 of the 
source condensate along direction i: 

~(i. \ - ^Po-R \2Ji(kjRj) - k i R i J (k i R i )\ 

P\ K i) — n 77. n \fi ' V 14 i 



{hR t f 



The ratio of w\ /w\ 



1.14 for the TF parabolic 
density profile (which can be contrasted with the larger 
value of a \ jo\ — \/2 in the case of a Gaussian 
profile; see below) is in good agreement with the re- 
sults of positive-P simulations of He* BEC collisions 
for relatively short collision times (< 25 /js) [11|, for 
which the obtained ratios ranged between 1.08 and 1.13 
depending on the correlation direction i. Similarly, 
good agreement is obtained when comparing the indi- 
vidual widths, in which case the positive-P results were 



,(CL) 



1.27w. 
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and 



..(BB) 



„(S) 



l),z 



u iS) 



w. 



(BB) 



1.18w 



(S) 



1 .39wy.z ■ The somewhat larger values of the 



y, ^correlation widths than the above analytic results are 
explained by the fact that the simulated source conden- 
sates were confined in the transverse (y and z) directions 
much more strongly than in the axial (x) direction, and 
therefore the actual GP ground-state profile (which was 
used as the initial condition) along y and z was interme- 
diate between a Gaussian and a TF parabola. 

The present analytic results and the positive-P results 
for 4 He* [ll| are generally in satisfactory agreement with 
the experimentally measured absolute widths [2| , except 
that the experimentally measured CL width was some- 
what smaller than the BB width, which is in contrast to 
the above analytic results and the positive-P results for 
4 He*. An obvious suspect for this discrepancy is the fact 
that the effective collision time in the experiment was 



~ 150 /us, whereas the positive-P simulations of Rcf. 11] 
were restricted to < 25 ps due to the large sampling 
errors developed past that time. Therefore, to explain 
this discrepancy it is important to perform first-principle 
simulations for longer collision durations and monitor the 
long-time dynamics of the correlation widths. 

To this end we have been able to perform such sim- 
ulations for collisions of 23 Na BECs, in which case 
the smaller scattering length and larger mass compared 
to 4 He* give more favorable parameter values for the 
positive-P simulations and allow us to extend them to 
collision durations < 650 /is. More specifically, we have 
performed the simulations as in Ref. [9(, starting from 
the full effective field-theory Hamiltonian, and extracted 
the CL and BB correlation widths as a function of time. 
The results are shown in Fig. [T] where the marked solid 
and dashed curves refer to the numerically obtained CL 
and BB widths, respectively. The horizontal solid and 
dashed lines, on the other hand, are the corresponding 
CL and BB widths from the present analytic treatment, 
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the same procedures as with the TF parabola, we obtain 
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FIG. 1: Momentum widths of the BB and CL correlations 
relative to the width of the source BEC in x, y, and z direc- 
tions as a function of time. Fhe solid and dashed curves with 
marks are the results of the positive-P simulations for 23 Na 
BECs as in Ref . Q . Fhe horizontal solid and dashed lines are 
the corresponding correlation widths from Eqs. (|12|l and (|13|) . 



Eqs. (fl2|) and (|T3|) . We see that the short-time asymp- 
totic limits of the exact numerical results converge to the 
analytic predictions, thus proving the validity and use- 
fulness of the developed analytic approach in this limit. 

The long time dynamics of the correlation widths from 
the positive-P results show that at some point in time 
the BB correlation width becomes larger than the CL 
width. This supports the experimental observations of 
Ref. [2J that Wy, z >w VtZ . It is interesting to note that 
the crossover point (~ 300 /is) is close to the point in 
time (~ 400 /is) where the peak BB correlation drops 
below the collinear HBT peak value of 2 (see Fig. 3(a) 
of [9(]). Both these observations seem to be unique to the 
exact first-principle simulations and have not been ob- 
served so far using approximate theoretical techniques. 
We attribute these effects to the dynamical broadening 
of the momentum distributions of the expanding conden- 
sates, and partly to rescattering of scattered atoms with 
opposite momenta back into the condensate modes. 

(21 

The importance of the observation of g^i(Q, i) < 

( 21 

<7cl(0>*) is related to the prospect of observing relative 
number squeezing between the s-wave scattered atoms 
with equal but opposite momenta 1 1] . The relative num- 
ber squeezing itself is a manifestation of a violation of 
the classical Cauchy-Schwartz inequality, corresponding 
to (?bb(0, t) > 3ql(0, t) [14J, while the opposite inequality 
corresponds to absence of squeezing. While we certainly 
do have a violation and relative number squeezing in the 
short-time limit 11], the above observation that #bb(0, t) 

(2) 

becomes smaller than g^(0, t) implies that the squeezing 
is lost for long collision durations. 

Gaussian density profile. -- We now give the results 
for a Gaussian density profile of the source conden- 
sate, ;o(x) = poGxp(— J2i x i/^ a i)' corresponding to the 
momentum distribution of p(k) oc exp(— ^. kf/2a k ), 
where o~i and a^ = l/2<7i are the rms widths. Following 



.,(2) 



g££(Aki,t « to) ~ 1 + exp(-Afc//8a^ 



(15) 



fl^(Afcj,t«*o) 



*V(o) 



exp(-Afc?/4tff). (16) 



r (CL) 



The correlation widths are now given by o~ k . = 2(7^ 
and a k . = V^c^, resulting in the ratio a k . /cr k . = 
\/2. These results are in agreement with those of Ref. [10( 
(and Ref. [6j under the same approximations as here) and 
represent an alternative derivation in the Gaussian case. 

In summary, we have developed a simple perturbative 
approach to pair correlations between the atoms scat- 
tered from a collision of two BECs. The analytic results 
obtained show how the CL and BB correlations depend 
on the shape and the size of the colliding condensates. 
The results are compared with exact positive-P simula- 
tions and are in good agreement in the short time limit. 
The long time dynamics - accessible through the positive- 
P simulations of sodium condensates - give evidence that 
the BB correlation width grows with time faster than the 
CL correlation width and that it can become broader 
than the CL width. This conclusion agrees with the ex- 
perimental observation of Ref. [2] , in contrast to previous 
theoretical predictions. 
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